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Abstract: The simplest extension of the super symmetric standard model - the addition 
of one singlet superfield - can have a profound impact on the Higgs and its decays. We 
perform a general operator analysis of this scenario, focusing on the phenomenologically 
distinct scenarios that can arise, and not restricting the scope to the narrow framework of 
the NMSSM. We reexamine decays to four b quarks and four r's, finding that they are still 
generally viable, but at the edge of LEP limits. We find a broad set of Higgs decay modes, 
some new, including those with four gluon final states, as well as more general six and eight 
parton final states. We find the phenomenology of these scenarios is dramatically impacted 
by operators typically ignored, specifically those arising from D-terms in the hidden sector, 
and those arising from weak-scale colored fields. In addition to sensitivity of mz, there are 
potential tunings of other aspects of the spectrum. In spite of this, these models can be very 
natural, with light stops and a Higgs as light as 82 GeV. These scenarios motivate further 
analyses of LEP data as well as studies of the detection capabilities of future colliders to the 
new decay channels presented. 
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1. Introduction 

Since 1934 when Fermi wrote down his theory of the weak interactions, we have singled out 
the weak scale as an interesting scale for new physics. Within the last two decades, we have 
finally reached this energy scale at LEP and the Tevatron, and with the LHC we should be 
able to probe this scale thoroughly. 

The most pressing question at the weak scale is the origin of electroweak symmetry 
breaking (EWSB). Within the Glashow-Salam- Weinberg theory of weak interactions, it is 
broken by a fundamental scalar doublet of SU(2), the Higgs field. In this framework, all 
precision quantities are calculable and agree with present experimental limits. 1 

The major theoretical shortcoming of the standard model (SM) is the question of the sta- 
bility of the Higgs mass. Scalars in general acquire quadratically divergent mass, suggesting 
that new physics should cut off this divergence and appears near the weak scale. Numer- 
ous solutions to this theoretical problem have been proposed, most notably supersymmetry 
(SUSY). 

We are of course neglecting dark matter here, neutrino masses, the LSND anomaly and the NuTeV 
anomaly, all of which are beyond the scope of this paper. 
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2. Naturalness in the MSSM 



In this paper we shall focus on supersymmetric solutions of the hierarchy problem. In these 
theories the introduction of superpartners cancels the quadratic divergences of the standard 
model with loops of opposite statistics particles, leaving only logarithmically divergent con- 
tributions to the Higgs mass, proportional to the SUSY breaking soft mass parameters. 

Despite its excellent success in controlling divergences from very high scales, in recent 
years SUSY has become far more constrained, with superpartners pushed to higher scales, 
resulting in tunings which are typically 0(4%) in the minimal supersymmetric extension of 
the standard model (MSSM). 2 The reason for this is quite simple: the success of the LEP 
program has pushed the lower limit on a SM Higgs boson above 114 GeV, a limit which 
applies in a large region of the MSSM parameter space. In the MSSM, the Higgs has a mass 
which has a tree level upper limit of mz = 91 GeV. To achieve a mass of 115 GeV, one 
must invoke large radiative corrections due to (s)top loops. At one loop, the correction to 
the Higgs mass is [2] 



requiring a stop mass of roughly 500 GeV. Unfortunately, such a large stop mass also feeds 
into the soft mass squared of the Higgs, 



where A is the scale where the stop masses are generated. Even for low scale mass generation 
such as gauge mediation, where A^167r 2 m^, this results in a radiative correction to the 
up-type Higgs soft mass squared of O((400 GeV) 2 ). This must then be cancelled against a 
positive mass squared, for instance via a /U-term, in order to achieve the appropriate vev, 
v = 174 GeV. This cancellation must be tuned at the level of 4% in order to achieve this. If 
the stops are heavier, or the mediation scale is higher, the tuning becomes worse. 

The root of this problem is the small quartic of supersymmetric models, which is fixed 
by the electroweak gauge couplings. There are many proposals to enlarge the quartics, for 
instance by adding terms in the superpotential [3-10], with non-decoupling D-term quartics 
[11-13], or through strong dynamics at an intermediate scale [14-16]. 

New Higgs Decays and Naturalness 

An alternative approach is to evade the LEP limits indirectly. For instance, if the Higgs field 
decays in a non-standard way, LEP may not have been as sensitive to the decay. This is not 
as trivial as it may appear, in particular because limits on invisible Higgs decays are nearly 
as strong as those of b-quarks. If these non-standard decays are present in nature, there must 
be new states lighter than the Higgs boson, which themselves decay into ordinary SM fields. 
The simplest possibility is the presence of an additional singlet superfield to which the Higgs 

2 For an excellent discussion of the details of fine tunings in SUSY models, see [f ] . 




v 2 yf sin 4 p log^mr /mf), 



(2.1) 



6m Hu = -^IK +r4)log(A/m f ), 



(2.2) 
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can decay, and which then decays into SM particles, such as b-jets. In this vein, it has been 
recently argued that within the NMSSM, where a singlet field S acquires a vev to supply the 
//-term of the Higgs sector, certain non-standard decays are possible, leading to more natural 
theories [17]. 

However, there is still a strong need for additional analyses, for many reasons. 

• New combined limits from LEP [18, 19] exclude most of the generic parameter space of 
these models. Typically, the decay h — > 2a — > 4b has been used to allow lighter Higgses. 
However, new LEP combined analyses have basically raised the minimum Higgs mass 
for this process to 110 GeV, nearly as strong as the limit on the standard model Higgs. 
This raises the question of whether other non-standard decays can occur generically, to 
which LEP analyses would have been less sensitive. We will discuss these limits further 
in the next section. 

• The NMSSM is not a fully general scenario. In the NMSSM, the singlet acquires a 
vacuum expectation value to generate a [i-term. One typically requires that this be 
the true vacuum of the theory, and additionally makes assumptions of the form of the 
theory in the ultraviolet (UV). All of these things can distract from - and are beside 
the point of - the basic phenomenological questions relating to Higgs decays. 

• The sensitivity of mz to UV parameters is not the only measure of naturalness in these 
theories. Although a significant reason to consider these decays is precisely to ameliorate 
this tuning, one often ends up with additional tunings in order to evade experimental 
limits. Thus, while mz may be relatively natural other elements of the spectrum may 
be highly tuned, and necessary, to evade experimental limits. 

• The parameter space of new operators has not been fully explored. In the MSSM, the 
full set of soft breaking operators is considered. However, in adding a singlet to the 
theory, one can consider both the effects of D-terms in the supersymmetry breaking 
sector, and non-degenerate soft masses for the scalar and pseudoscalar components of 
the singlet, both of which can have significant phenomenological consequences. 

• The sensitivity of the decays of the pseudoscalar to the presence of new fields has 
essentially been ignored within the context of supersymmetric theories. Such fields can 
naturally induce the decay of the pseudoscalar to two glue jets, which have far weaker 
constraints. 

There are two important additional points to be made with regard to naturalness. The 
first is: what do we mean by tuning in the scenarios we are considering? Since we are interested 
in the low-energy phenomenological theory, it is impossible to quantify the sensitivity of mz 
to the UV parameters. However, we know that the tuning of mz arises in general due to the 
large values of the stop squark masses. Therefore, we shall use the stop masses as a proxy for 
this tuning. However, there is typically a tuning necessary to achieve the proper spectrum for 
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the lighter scalar states, and this is usually the most severe tuning in the models. Therefore, 
we shall quote both stop masses, as well as the scalar mass spectrum tuning for every model 
considered. 

The second point is: how natural is it to include a soft supersymmetry breaking operator, 
without an associated superpotential operator? In particular, how reasonable is it to include 
a trilinear scalar potential A-term operator without associated Yukawas? This question has 
been studied previously in the context of sterile neutrino masses in supersymmetry [20,21]. 
In fact, if there is a SUSY breaking field which carries a charge, whether R-charge or Peccei- 
Quinn charge, it is quite natural for these operators to appear without an associated superpo- 
tential term. The converse, however, is not true. The presence of the superpotential term will 
radiatively generate the soft term at a minimum at the loop-suppressed level. Consequently, 
we will insist on technical naturalness, that these soft operators are present at least at this 
small level. 

In lieu of these points, we will pursue a phenomenological study of the effects of singlets 
on the decays and properties of Higgs bosons. We will focus on decay modes that have not 
been considered previously, including final states with 6 and 8 particles. We will often be 
studying situations when both the new scalar and pseudoscalar states are lighter than the 
Higgs, and the Higgs is light enough to have been produced through Higgs-strahlung at LEP, 
although we will consider other scenarios. By also focusing on the relevant parameters, we 
hope to clearly elucidate the effects of mixing of the singlet with the Higgs, and the effect on 
the mass of the Higgs boson. 

Since this paper has both basic phenomenological points, as well as points related to model 
building, we have attempted to lay out this paper so that one who is interested principally 
in the phenomenology can still learn the relevant points. Sections headed "Model Building" 
are independent, and the paper can largely be read without those. 

The layout of the paper is as follows: in the next section, we will review existing limits on 
Higgs decays, both for SM-like and non-SM-like Higgs bosons. In section ||, we will discuss the 
effects of singlets on Higgs physics. First, we show how non-standard decays into scalars can 
arise and dominate the Higgs decay width. We will also point out how non-standard decays 
allow a larger mixing with the Higgs, this mixing can raise the mass of the Higgs boson 
considerably without resorting to radiative corrections. In its "Model Building" subsection, 
we present the relevant operators that induce the Higgs mixing and decays. In section [|, we 
will sketch out the spectra and non-standard decay modes which are consistent with existing 
LEP bounds, and lead to a more natural parameter space of the theory. In its "Model 
Building" subsections, we discuss model realizations of these different spectra and decay 
scenarios. In doing so, we can assess more simply the degree of tuning required to satisfy 
the experimental constraints. In section |(| we summarize the presented benchmark points 
and their relevant phenomenology for future collider experiments. In section |^, we discuss 
additional directions worthy of further investigation and conclude. In particular, we suggest 
some new analyses on LEP data which might be useful when considering these scenarios. 
Finally, in appendix [A] some calculational results for scalar trilinears are summarized. 
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3. Summary of LEP limits on Higgs 

Most particle physicists are familiar with the LEP2 95% CL limit on the SM Higgs of rrih > 
114.4 GeV, but LEP has also produced a multitude of other limits on Higgs physics. 

The limits that are most applicable to this work are the 95% CL limits on the so called 
£ 2 parameters (often also referred to as C 2 ). We are analyzing a SM-like Higgs, which means 
that Higgs-strahlung production of the Higgs is close to the SM rate. 3 Then depending on the 
assumed process of Higgs decay h — > X, the defined limits are on £ 2 = .(e+ e -->h > z^M -^-^(^ — * 
X). We now list the different limits that LEP has analyzed, based on Higgs-strahlung and 
the given Higgs decay process (note: all stated mass limits assume £ 2 = 1): 

• Model Independent Decays: This is the most conservative limit on the Higgs boson. It 
assumes that the Higgs is produced with a Z boson and looks for electrons and muons 
that reconstruct to a Z mass, while the Higgs decay process is unconstrained (by theory 
or the event analysis). The only study of this sort is done by OPAL giving a limit of 
rrih > 82 GeV, see Fig. 11 of [22]. Unfortunately, no other collaboration has released 
such an analysis. 

• Standard Model Higgs: LEP-wide limits on the SM Higgs are given in Fig. 10 of [23], 
requiring rrih > 114.4 GeV. This study also includes the strongest limits on h — ► bb,rf 
rates. 

• Invisible Decays: In this analysis, the Higgs is assumed to decay into stable (on collider 
length scales) neutral particles. The implication is that non-standard Higgs decays have 
to primarily decay into visible particles. Both L3 and DELPHI have performed such 
an analysis [24, 25] , but the most stringent constraints are from an older preliminary 
LEP-wide analysis that has a limit rrih > 114 GeV, see Fig. 4 in [26]. 

• 2 Photon Decays (aka Fermiophobic) : Fermiophobic typically means a Higgs with stan- 
dard couplings to gauge bosons, but suppressed couplings to fermions, allowing decays 
into WW*, ZZ* as well as photons. If all decay modes are open, there is a limit of 
mh > 109.7 GeV, while decays exclusively to two photons have a limit of rrih > 117 GeV. 
See Fig. 2 of [27], which is the LEP-wide analysis. 

• 2 parton hadronic states (aka Flavor-independent): In this type of analysis, any two jet 
decays of the Higgs are allowed. The analyses use the jets that are least sensitive to the 
candidate Higgs mass and details of the Z decay. Each LEP collaboration has done a 
study [28-31]. However, the strongest constraint is the preliminary LEP-wide analysis 
of m h > 113 GeV, see Fig. 2 of [32]. 

3 More generally, we will usually be taking the decoupling limit, where the lightest Higgs state is produced by 
Higgs-strahlung, but not in associated production with the CP-odd A . Thus, constraints from Higgs-strahlung 
are the main concern in these scenarios. 
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• Cascade Decays: These are the constraints that are most relevant for the present study, 
where cascade decays mean that the Higgs decays into two scalars <f> and those scalars 
decay into X (i.e. h -» 2<j) -> 2X). OPAL [33] and DELPHI [34] looked at b decays 
(X = bb), see their Fig. 12's and a new LEP-wide analysis [18,19] has constrained both 
b and r decays, with the exclusion plots given on page 8 in the first reference. For 
h — ► 2(f) — ► 46 the limits are now 110 GeV for a Higgs produced with SM strength. For 
other intermediate scalar decays, 4> ~~ * 2g, cc, rf, the best model independent exclusions 
are from OPAL's analysis when the mass of the scalar is below bb threshold, which is 
given in Fig. 7 of [35] (note: the analysis is restricted for Higgs masses in the range 
45 - 86 GeV). 

4. Singlets and the MSSM 

As we have previously said, the simplest extension of the MSSM is to add a SM singlet 
superfield S, containing a CP-even scalar s, and a CP-odd pseudoscalar a. 4 The dominant 
phenomenological effects of this new field are: 

• New decays for the Higgs boson. If one or both of the new states are lighter than half 
the Higgs mass, decays h — » 2s and h — ► 2a are possible, followed by s — > X or a — > X. 

• Light states can mix with the Higgs boson. If a light state mixes with the Higgs boson, 
it can push the mass of the Higgs boson up without large radiative corrections. This is 
at the cost of having a new light state which can be produced through brehmstrahlung 
off a Z, for which there are stringent constraints. 

We first show that non-standard decays can easily dominate over SM decays. In the case 
of cascade decays, for a Higgs produced with SM strength through Higgs-strahlung, the new 
process must have a width at least 5 F/ l _ >2 a ~ 4 x r^ 2 6- With a term in the Lagrangian 

^=ha 2 (4.1) 



V2 



c 2 v 2 ( Am 2 ^ 



the width to two scalars is [36] 

T h ^ 2a = — — — ( 1 - ) (4.2) 
167rm fc V m h J 

which one can compare with the dominant width of the Higgs to b quarks (in this mass range) 

N c ml f 4ml\ 3/2 nnn „^ T ( m h 



4 More precisely, the CP-even properties of a are not determined until its interactions and mixings are 
given. For example, it can pick up a CP transformation if it mixes with a scalar with fixed CP properties (the 
CP-even h° , H° or CP-odd A ) or couples to fermions/gauge bosons in a certain CP fashion. On the other 
hand, in all cases s mixes with h° and thus is always CP-even. 

5 This condition is not sufficient if the scalar a decays into b quarks, as will be discussed later in section 5.f 
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Taking into account higher order effects, a more accurate approximation for the Higgs decay 
width to SM particles in the mass range of interest is (using HDECAY [37]) 

Thr-tSM = (2-34 x W~ 5 m h + .206 MeV) (4.4) 

For a 100 GeV Higgs, in order for scalar decays to sufficiently dominate, one must have 

^= > 5 GeV (4.5) 



The required size of the trilinear indicates that non-standard cascade decays can easily dom- 
inate over standard Higgs decays. 

One can also study the effects of mixing with a singlet on the Higgs mass in a model 
independent fashion. In the presence of mixing, the mass eigenstates s and h will be related 
to the interaction eigenstates through a mixing matrix 



s \ ( cos 6 — sin 6\ ( s 
h ) V sin# cos# j \ h 



(4.6) 



Similarly for the mixing in the CP-odd sector, the mass eigenstates a and A are related to 
the interaction eigenstates by 

'a\ = fcos(j) -sin<A / a\ 
A ) \sm$ cos^ / / 

If the MSSM mass of the Higgs is , m'i nssrn , then the mass eigenvalue after mixing with the 
lighter singlet is 

m\ = m ™~ 9 m J Sin2g , (4.8) 
" cos^ 9 

where mg is the mass of the mass eigenstate s. This increase in mass is only through mixing 
and not through radiative corrections. 

In the subsection below, we discuss some of the operators which can be used in building 
models with certain phenomenologies. Those readers who are only interested in the descrip- 
tions of the phenomenology can proceed to section H. 



Model Building: New Operators with Singlets 

The introduction of a singlet superfield allows us to introduce a number of new operators in 
the theory. These operators have the effects of introducing into the theory masses, mixing, 
and decays of the states. The operators are listed in Table |l[ We will describe here the 
effects of these operators when they appear singly. Clearly, the effects of the operators do not 
necessarily add linearly, particularly when vevs turn on, but it is useful to understand their 
effects individually. In many scenarios this is sufficient to understand the phenomenology and 
we will point out the situations where there are interference effects from multiple operators. 
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6 2 s 


XXt(S + S*) 2 D 


mass for s 


8 2 


xxUs - S^n 


mass for ci 


m D 


{W' a W a S) F 


singlino (Dirac) mass, mixes s and h 
allows h — > Is decays, allows Z s production 




{SH u H d ) F 


(with mu term) mixes s and h, s and a masses, 
singlino (Dirac) mass, allows h — ► 2s, h — ► 2o decays 


A 


jp + tl.C. 


Real: (with (s)) mass for s s — > 2a decays, 3 s coupling 

I in a cri n a • // ^ 9c Hpipsnrt; s / / fnnn huff 

LLLLCX/tl LLLCXiL y . (. ( ? £jO ucua V a , O (/ K^yj Ll IJllilli 


A h 


(XSH u H d ) F + h.c. 


mixes s with H, a with ^4°, a decays (due to mixing) 
with mixing, allows h — > 2s, 2a decays (esp. large tan/3) 


m CP 


XXi(S + &)(S-&) D 


Mixes s and a, allows a decays (due to mixing) 


M Q 


(SQQ + MqQQ) f 


Allows s and a decays to gluons and photons 



Table 1: Operators involving the Singlet 



S 2 , S 2 — Scalar and Pseudoscalar Soft Masses 

The operators arising from Kahler potential terms, S 2 and 5 2 , determine the spectrum of the 
theory by adding mass terms to the Lagrangian 

£ 5s2 = -|s 2 Ag = -|« 2 . (4.9) 

Most previous analyses include equal soft masses for the scalar and pseudoscalar, although 
there is often no a priori symmetry reason to do so. All other scalar fields in the MSSM 
carry a gauged U(l), forcing a degeneracy which is not necessary for singlets. Because our 
interests are phenomenological, we allow different masses, and thus more interesting spectra. 
In particular, this allows spectra of the form rrih > 2m s > 4m a , leading to events with many 
final state jets, b-quarks and photons, produced through cascades. 

The tuning of these theories with singlets, when < 114 GeV, is usually encoded in 
how tuned these additional masses must be in order to achieve a proper spectrum, and not 
in the sensitivity of mz or the masses of the stops. When we refer to these operators, a value 
of 5 2 or 5 2 of ±(20GeV) 2 refers to a contribution of ±(20GeV) 2 to the scalar or pseudoscalar 
mass squared. 

mo - Supersoft Operator 

m.D is the so-called "supersoft" operator [38] and arises from a SUSY breaking spurion in the 
superpotential. This operator is particularly interesting, and often neglected. If the hidden 
sector has a U(l)' which acquires a D-term expectation value, it has no effect in the MSSM, 
in that all soft operators can be generated by F-term breaking. On the other hand, in the 
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presence of a singlet, D-term breaking becomes important through the introduction of this 
new operator [39]. The Lagrangian terms generated are 



f o W 

= / d 2 9^W?S + h.c. 



M 



-(itids + D Y y 



(4.10) 



where Dy is the usual hypercharge D-term, Dy = '}2 i gyqi<j)*4>i, and -jf = 6 a rriD- Although 
it is supersymmetry breaking, it does not feed into the RG flow of any other soft operators, 
hence the term supersoft. The principal effects are to produce a mass for s (but not a), and a 
coupling shh* which is proportional to the new contribution to the s mass. It also generates a 
Dirac mass for the Bino with ips equal to m£>/2. When we refer to a value of mo of 20 GeV, 
this implies a contribution to the scalar mass squared of (20GeV) 2 , and the corresponding 
trilinear term. 

This operator directly affects the physics of electroweak symmetry breaking and the 
properties of the Higgs. However, it is simple to perform the same minimization of the Higgs 
potential in the presence of this term. Let us consider the usual MSSM neutral Higgs potential 
with this contribution. 



V = (M 2 + m 2 Hu )H 2 u + (H 2 + m 2 Hd )H 2 d - {bH u H d + h.c.) 

r,1 



8 
1 

+ 2 



(4.11) 



m DS + ^(H* u H u -H* d H d ) 



A 2 

+ T s2 



where A 2 encodes other contributions to the s mass, such as from 5 2 . 

When electroweak symmetry is broken, there is a linear term for s, and so we shift s by 
an amount 

gym D v 2 cos 2/3 

The Higgs potential (neglecting s which has already been shifted to its appropriate minimum) 
now reads 



V = 



H 2 + m 2 Hu + 



g Y m 2 D v 2 cos 2j3 



4(m 2 D + A 2 ) J u 



M 2 + ™L" 



g Y m 2 D v 2 cos 2/3 
4(m 2 D + A 2 ) . 



H 2 (4.13) 



(bH u H d + h.c.) + + / y) {H* U H U - H* d H d ) 2 



So we see the presence of the s vev amounts to a redefinition of the values of rn 2 Hu and m 2 Hd . 
Thus, we can perform the usual minimization and diagonalization of the MSSM Higgs fields. 
The new trilinear couplings induce a mixing between the s and the other Higgses, which is 
encoded in the mass matrix (in the (s h H) basis) 



m 2 D + A 2 s 

9Ym D vs a+/3 

gYm D vc a+fj 
V2 



gY'rriDVs a+l3 


g Y m D vc a+[) 






m\ 








m H 



(4.14) 
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Here, as we dial up m 2 D , which is always positive, we increase the mixing. The values of m\ 
and m 2 H are the usual ones from the MSSM, in particular with m\ bounded at tree level to 
be below m|. 

Since s a+ p = sin(a + 0) > cos(a + f3) = c a+ p generally, and m 2 H > m\, we can just focus 
on the 2x2 mixing submatrix between s and h. The mass eigenstates are related to the 



interaction eigenstates as given in (4.6). We are principally interested in the trilinear terms, 
especially the s 2 h term, which allows h —* 2s decays. There are two contributions to this 
term, one arising from the new supersoft operator, once we have mixed the CP-even states, 
and one from the usual D-term trilinears, once we have mixed. The coefficients of these terms 
are given in the appendix. 

A^SHuHd 

Xh is a commonly studied superpotential operator, as it induces an additional Higgs quartic 
at small tan/3 which can raise the tree level Higgs mass [8-10]. After electroweak symmetry 
breaking, it generates masses for a and s of size XhV. With these contributions the potential 
is given by 



V = (H 2 + m 2 Hu )H 2 u + (H 2 + m 2 Hd )H 2 - (bH u H d + h.c.) + 
i9 2 +9 2 y), 



■{H* U H U - H* d H d y + X z h (H u H d )(H u H d y (4.15) 



-^(s 2 + a 2 )(H* u H u + H* d H d ) + V2X hf is(H* u H u + H* d H d ) + |L 2 + |l 



Note the explicit //-term for the Higgses (and the resulting trilinears), this is to be compared 
to the NMSSM where the vev of S gives the //-term. However, absent any additional mass 
terms for S, i.e. 5 S , we can simply shift S — > S — y/2[i/Xh, removing the entire \x term. In 
the presence of of additional soft masses for s there is an s-vev of, 

<•> = -=T$4- <« 6 > 

The effect of this is to replace /i by an effective n, which we denote jl, where 

Through the new quartic and trilinear terms there are mixings between s and the MSSM 
Higgses, giving a mass matrix 

^W + ^s -2X h vfis a _f3 2X h vflc a ^p\ 
-2X h v{is a _p m 2 h (4.18) 

2X h vjlc a -.p m\ j 

Notice, however, that in the absence of additional soft masses for s (5 2 = 0) the mixing 
vanishes, and increases as we deviate from zero. The presence of s-h mixing allows s — > 2a 
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decays. It is possible for h — > 2s, 2a decays to occur with equal amplitude, for instance at 
small mixing. However, it is also possible to arrange for A/j to give only one dominate decay 
of h. Finally, remember that the presence of A^ naturally requires the presence of Ah, at least 
at a loop-suppressed level. 

A S S S + c.c. 

Another possible operator to add is a soft trilinear for S in the potential, AgS 3 + c.c. Alone 
this has no effect on the Higgs, although it does give equal but opposite contributions to s 
and a soft masses with an S vev and also induces s — > 2a decays. However, if there is another 
source of s — h mixing, A s can generate h — > 2s and h — ► 2a decays as well as the original 
s — ► 2a. The effective potential with this contribution is, 

V = (H 2 + m 2 Hu )H 2 u + (H 2 + m 2 Hd )Hj - (bH u H d + h.c.) + (4.19) 

{ -l±tt±(H:H u - H* d H d f + ^ - 3,a 2 ) + |, 2 + fa 2 

Since this operator does little of interest by itself it must be considered in tandem with 
some other operator that produces s — h mixing and an s-vev. One combination of note is 
mo and A s . We can analyze the general case under the simplifying assumption, as above, 
that the mixing involves mainly s and h and not H, so the mass and interaction bases are 
related as in Eq. (|4.6|) . 

The relevant trilinears are given in the appendix. As in the case of A/j, at small mixing 
h — ► 2s and h — > 2a have equal amplitudes. The size of the s-vev (s), generated by the mixing 
operator will be corrected by A s but if A s (s) 2 <C m? s (s), as will often be the case, its effects 
can be ignored. 

AhSHuHd + c.c. 

The potential with Ah included is, 

V = (M 2 + m 2 Hu )H 2 u + (H 2 + m 2 Hd )H 2 - (bH u H d + h.c.) + (4.20) 



(9 2 + 9y) {KHu _ H * Hd) 2 + {AhSHuHd + c . c) + 6 1 S 2 + 6 J a 2 



2 2 

After EWSB this generates a vev for s, 

Auv 2 

(°) = -jt r V (4.21) 

This gives a correction to the b-term for the MSSM Higgses, but still allows the usual diago- 
nalization, leading to a mixing matrix of the form, 

S 2 A h vc a+/3 A h vs a+P s 
A h vc a+P m 2 h ) ( 1.22) 

, A h vs a+I3 m 2 H 
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Unlike the previously discussed operators, Ah mixes the CP odd component of S with A . 
The mixing matrix, in the basis (a A ), is 



(4.23) 



1117, ~Ahv\ 



-A h v m\ J 



The mass eigenstates are related to the interaction eigenstates as given in ( 4.27| ). This mixing 



of a with ^4° means that a can decay to 2 b-quarks or 2 r's if its mass is below about 10 GeV. 
With no additional terms in the Lagrangian, e.g. Mq 1 ,rr? CP = 0, this is the only decay path 
of a. The relevant trilinears are given in the appendix. 

rr? CP s a — CP Mixing Mass 

The operator CP mixing mass which mixes the would-be pseudoscalar a with the 

would-be scalar s. This operator can be important in allowing decays of the a when Mq 
and Ah are zero and can be particularly important when decays h — > 2s would dominate the 
Higgs decay, but are kinematically forbidden. In this case, mixing due to rr? cp can make 
h — > sa the dominant decay mode. 

We should point out that this operator does not, by itself, violate CP. As mentioned 
earlier, without any other interactions, there is nothing to prevent us from assigning CP-even 
transformation properties to both a and s. This remains true with the inclusion of A s or the 
supersoft operator wid in the potential. Only when a couples to fermions or gauge bosons as a 
pseudoscalar, or mixes with the A is a CP-odd property forced upon it. The scenarios which 
we will study will not have such a mixing and hence no actual CP violation is introduced into 
the theory, making such a term safe from edm searches, for example. 

One requires either rri£, or Xh to induce Higgs decays (through direct couplings or mix- 
ings). In the case that A^ is present, h — > 2a is already significant, so h — > as typically cannot 
compete. However, in the mp case, where a has no couplings, such a scenario is viable. We 
thus consider the mass matrix (in the (a s h)) basis) 

$1 ™ 2 cp 




m 2 D + 6 2 9 ym D vs a+fj ^ 



m CP 

„ gYiri D vs a+l3 

Because the mixing will allow the a to be produced through o-strahlung from the Z, LEP 
limits become very severe. Ultimately, this implies that vr? CP must be somewhat small, and 
for our purposes here, we will treat it as a perturbation. Ultimately, it will be necessary 
to calculate mixings precisely, but for estimates and intuition, it is easiest to study in the 
perturbative limit. We can perform a rotation in the sh sector by an angle 8, which may be 
large. This leaves the following matrix 

m a m CP cos ^ m CP sm ' 

m 2 CP cos6 mj j (4.25) 

m~ 

h 
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The mass matrix can then be diagonalized by a series of rotations 



U = R(6 sh )R(5 as )R(5 ah ) (4.26) 

where 6 s h = diagonalizes the sh mass matrix, and the R(5i) are perturbative rotations 
where sin<5 as = m 2 CP cos /{ml — m\) and sin 5 a h = rricp sin0/(m~ — m 2 ), where m 2 ^ are the 
ss and hh entries after the 9 rotation. 

The relevant mixing angle which enters into a-strahlung is 

sin d a h = sin 9 cos 5 a h sin 5 as — cos 8 sin 5 a h — mcp sm ^ cos ^ ~~ m a)~ 1 ~ ( m | ~~ m a) _1 ) ■ 

(4.27) 

Mq 1 — Fermiophobic Decay Operator 

If the singlet superfield couples to new vector matter, there are loop induced decays of the 
scalars into gauge bosons. The presence of the new fields corrects the beta function of the 
SM gauge groups, with a superpotential coupling 

XqSQQ + M Q QQ (4.28) 

the gauge kinetic term becomes 

/ dH (2^ + 16^ bg(mQ + XqS) ) WaW ° " (2^ + -^ M Q S ) WaW ° (429) 

where W a is the gauge field strength, and Mq 1 = Am^ 1 . This results in a coupling between 
s and two gauge fields, as well as a and two gauge fields which explains why we refer to this 
as a fermiophobic decay operator. Because s generally mixes with the Higgs, it is difficult for 
a loop suppressed decay to compete. However, a, which mixes with the often heavy A , and 
often through the loop suppressed operator A^, can have its dominant decay mode through 
this operator. Expanding out the expression above, we find a term in the Lagrangian 

C D c(aF a F a ) = c{ae^P a F^F^)/2 (4.30) 

This induces a decays into photons or gluons, with the decay width 

T c = ^m 3 a N D (4.31) 

where Np is the number of "colors" in the final state (i.e. 1 for photons and 8 for gluons). 
With the convention above, 

c = ~64^ b ^ = ~[^r hai (432) 
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where i labels the gauge group and hi is it's beta function where d f^^ 
specific case of a 5 + 5 and higgsinos 



16tt 2 





D,D C 


L,L C 


H u , H d 


^3 


1 








bQ 


2/3 


2 


4/3 



In the 



(4.33) 



We can compare the decays into gluons and those into b quarks through mixing, 

(iV c 2 - l)c 2 ml/(4n) 



glue/^bb 



(nib I v) 2 iV c (sin <f> tan j3) 2 ^jl — Am 2 ./ 'm 2 a [m a / 'IQit) 
which for the particular case of a D,D C , becomes 



(4.34) 



glue/^bb 



4tt 2 



N c (sin 4> tan 0) 2 yjl — Am\jm 2 a 



m a v V 
2 \m b MQj 



(4.35) 



For decays to glue to dominate over b's, one requires small mixing angles in the CP-odd sector, 
which can be achieved simply with loop suppressed Ah, or a heavy A . Alternatively, one 
can have lighter D,D C squarks, for example, with the values 03 = .12, m a = 40GeV,MQ = 
330GeV one requires (sin^tan/3) ^$ .05 to get comparable rates. 

When Mq 1 induced decays for a dominate, we can get an estimate on the branching 
ratios into gluons and photons in the case of S coupling to a SU (5) complete multiplet. For 
this, the branching ratios are (2*7,27) = (.99,3.8 x 1CT 3 ). If the Higgs dominantly decays 
into 2a, this gives branching ratios h — > 2d — > (4g, 2g 27, 47) of (.99, 7.6 x 1CT 3 , 1.5 x 1CT 5 ). 
The channels with photons may be enhanced by including which only induces decays into 
photons, via Higgsino loops. As a final comment on this operator, we point out that if this is 
the only allowed a decay, as Mq increases, a will first decay with a visible displaced vertex and 
eventually will decay outside of the detector. Both of these possibilities are highly constrained 
by LEP searches if the Higgs can be produced, ruling out such a Higgs that cascade decays 
into two such a scalars. 6 Therefore, we only consider values of this operator where a decays 
promptly. 

Other Operators 

Note that we ignore here the effects of a /i-term for the S field, and a superpotential S 3 term. 
Since we are mostly concerned with the scalar sector of the theory, and aren't concerning 
ourselves with the cosmological implications at this time, the mass of the singlino is irrelevant 
for our purposes. The only effect of u, s is then to generate an S*H u Hd term in the potential 
in a cross term with A^. This has the same phenomenological effect as Ah, so including /is 



6 Note added: It has recently been emphasized that similar Higgs decays with highly displaced vertices 
could well be visible at Tevatron and LHC, enhancing its detection prospects [40] . 
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should not change the basic outcome of our study. Including S 3 in the superpotential can 
induce invisible decays of the Higgs to a singlino, but these are highly constrained and not 
interesting for our purposes. It can also induce a term H u HdS* 2 in the potential, but this 
acts simply like Ah in the presence of an s-vev, and like Ah in the presence of an h-vev in 
allowing h — > 2a, 2s decays. Hence, we do not consider either operator here. 

5. Scenarios 

Although we have expanded the MSSM only by a single superfield, there is a remarkably large 
number of new scenarios of Higgs decays which arise, with large variations in the number of 
particles in the final state, as well as types of particles in the final state. 

We can group the scenarios, in increasing order of complexity, into three categories: 

• The Higgs is mixed heavy. In the presence of large mixing, the Higgs mass can be 
increased significantly, in spite of small radiative corrections. In this scenario, the 
Higgs is at least 114 GeV in the case of SM-like decays, or 106 GeV in the case of 
cascade decays to 4 b quarks. 

• The Higgs decay is dominated by a single stage decay. Here h — ► 2a — > 2X or h — > 2s — > 
2X, where X is composed of a pair of standard model fields. As above, if X = 2b, one 
requires rrih 106GeV, however, if X = 2r or X = 2g, the Higgs can be considerably 
lighter. 

• The Higgs decay is dominated by a two-stage cascade. That is h — > 2s — > 4a — > 4X or 
h — ► as — > 3a — > 3X. Such two stage cascades generally do not occur if one restricts 
oneself to the scenario of the NMSSM. Here X = 2r, 26, 2g, resulting in 6 and 8 particle 
final states, which have not been constrained by LEP analyses other than the model 
independent ones. 

Additionally, we can split the scenarios further into two cases: "large-mixing" and "small- 
mixing". We define these as follows: the small mixing case arises when the light CP-even 
singlet has a sufficiently small Higgs component that s-strahlung limits do not constrain it, 
even with conventional decays. In the large-mixing case, the singlet has sufficient Higgs 
component that non-standard decays are needed to evade limits. 

With large mixing case, there are three real scalars that we are concerned with: h, s, 
and a. These are defined to be the field which couples to the Z in Higgs-strahlung, the field 
which mixes with the Higgs, and the field which does not mix with the Higgs, respectively. In 
many cases s can be thought of as the scalar and a as the pseudoscalar. The mass eigenstates 
are h, s and a, which are the fields mostly made of the untilded field of the same character, 
and usually the heaviest, intermediate, and lightest mass eigenstates, respectively. As a note 
of caution, we will often refer to the h, s, a still as the Higgs, CP-even singlet and CP-odd 
singlet and have tried to ensure that the meaning is clear given the context. 

First, we will discuss the necessary operators in all scenarios and then we will proceed to 
study the three basic scenarios, including the existing limits and model building possibilities. 
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Model Building: Necessary Operators in All Scenarios 

The basic requirements on any viable scenario are threefold, and can be satisfied with various 
combinations of operators. The requirements [operators satisfying the requirements] are: 

• a must decay [Ah, rn^p, Mq 1 ]: This is not absolutely true, for instance if a is never 
produced, but almost always is. 

• h must decay into 2s or 2d or as [mo, A/i, A s (+ mixing), Ah]: For the mass range we 
are interested in, these non-standard decays for h must exist. 

• s must decay to 2d [A^ (+ mixing), A s , Ah (+ mixing), m^ p ]: When s is light (i.e. 
i$ 80 GeV) and in "large-mixing" , constraints from SM decays are usually quite strin- 
gent. Conceivably, s could decay directly primarily to r's if the mixing is small enough, 
but in general, the cascade decays for s are necessary to evade LEP limits, unless it is 
heavier than h. 

5.1 Higgs Decays Through a Single Stage Cascade Into b Quarks 

As we have already described, the only means in the MSSM to push up the mass of the Higgs 
above LEP limits is to introduce very heavy top squarks, which then introduce tunings into 
the theory. However, as described in section f|, if the Higgs mixes with a lighter field s, the 
mass of the field h which is produced strongly through /i-strahlung from a Z can have its 
mass pushed up, simply by mass mixing, without introducing unnaturally heavy top squarks. 
However, this mixing allows the lighter field to be produced weakly through s-strahlung, 
trading off the heavier h for stronger limits on the lighter s. 

Now that s contains a Higgs component, we must consider LEP limits on it. With 
standard model decays for s, the limits are typically sin 2 9 ^ 0.03. The large mixing necessary 
to achieve changes in the h mass requires that s decays dominantly through cascades, in this 
case, one typically requires sin 2 9^,0.2 if the cascades end in b quarks for lighter s and 
sin 2 9 ^5 0.4 for heavier s. These, in general, also imply that the dominant decays of the Higgs 
are non-standard. 

If these decays are h — > At or h —* 4g, then the Higgs can be considerably lighter than 
115 GeV, and it is essentially unnecessary to push the Higgs mass up. These scenarios will 



be discussed in section 5.2 



We show the allowed regions (dark and light shaded) with h decays to four b jets in 
figures ||a and |l|b. In these plots, we assume a value for m a > 12GeV (where kinematics no 
longer strongly favor decays to r's) where constraints are strongest. We show the regions 
which are allowed when accounting for constraints on s — > 46 decays (light and dark shaded 
regions), assuming a typical BR a ^2b — 0.87. For comparison, we also plot the region where 
rrih > 110 GeV, where h — > 46 is essentially unconstrained (dark shaded regions). 

These plots have assumed a 100% branching ratio for h — > 2a, so the only applicable 
limits are the h — > 46 exclusions given in the LEP-wide analysis [18,19]. However, it is an 
important question what the LEP constraints are on a scenario with both h — ► 26 and h — > 46 
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Figure 1: The allowed regions in sin 2 9-mg space, where m$ is the mass of a light singlet mixing 
with the Higgs an amount sin#. The plots are shown with ^/m^ m ( - 2 = 270 GeV (a) and ^/mflTT^ = 
300 GeV (b). The shaded regions (combined dark and light) are allowed from the combined LEP 
limits, while the dark shaded regions additionally have > 110 GeV. 



rates. In a recent paper on these decays in the NMSSM [41], the authors only required that 
the h — ► 26 (46) rate be consistent with the LEP h — > 26 (46) exclusions, thus they assumed 
that these analyses are independent of each other. However, upon closer inspection of the 
details of the analyses, it appears such an approach may be too generous. 

In the higher mass (m; > 90GeV) range, the combined LEP analyses are dominated by 
OPAL [33] and DELPHI [34]. A close reading of these papers suggests that the 26 and 46 
decays of the Higgs are generally reconstructed together. For instance, in studying e + e~ — > 
Zh, in the Z — > vv analysis by OPAL, the neural nets trained to capture to the 46 decays 
of the Higgs are also reasonably efficient in capturing the 26 decays. In the all jet [Z — > jj) 
analysis, OPAL uses the same analysis procedure for both 26 and 46, forcing the six jet 
event into a four jet topology via the DURHAM algorithm, hence this analysis efficiently 
reconstructs both types of decays. As for DELPHI, they clearly state that the analyses are 
not independent of each other [42], and both decays are reconstructed via the same analysis 
procedure. 

Beyond efficiencies, there are still differences between 26 and 46 events. For example, the 
distribution of discriminating variables will be broader in 46 decays (e.g. the reconstructed 
Higgs mass), making it difficult to know how to constrain the scenario when there are signif- 
icant levels of both 26 and 46 decays. If 26 and 46 signals were indistinguishable, the correct 
limit would be — Y ^ 2b < 1, where £r d is the experimental bound of the individual anal- 
ysis. As an attempt to combine these limits, when there are O(l) rates for 26 and 46 decays, 
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in addition to applying the individual limits, we will also require that 



e 2b+4b ^^ + ^<V2 (5.i) 

$4b,bd $2b,bd 

This additional requirement, on the effective £| 6+46 , accounts for the redundancy of the analy- 
ses, and which should forbid situations where the naive combination of analyses is significantly 
excluded. Thus, we find it to be a good compromise between the assumptions of complete 
independence/interdependence (£26+46 — 2, 1 respectively) of the separate analyses. Note 
that one can also interpret this as only requiring a 99.5% CL limit, if the two decays give 
indistinguishable events. More rigorously, a combined analysis should be done to find the 
proper constraints. 

For our purposes here, where we take the branching ratio h — > 2a to be one, the results 
are relatively simple. Even with relatively light stops (270 GeV), the small suppressions of the 
Higgs couplings due to mixing are sufficient to allow such a Higgs to have been undetected. 
However, at this mass, this is achieved by strongly mixing and pushing the physical Higgs 
above 110 GeV(out of the LEP constraints). At mi = 300 GeV, a region opens up below 110 
GeV, but the overall allowed space is still quite narrow. Moving to higher Higgs mass (roughly 
325 GeV stops) the model independent parameter space opens up significantly (about twice 
as large as at 300 GeV). However, even with 300 GeV stops, the tuning of mz is already 
expected to be 0(15%), which begins to reintroduce fine tuning from another direction. 

Note that in all of these cases we are considering s to be lighter than h. The reason for 
this is simple: although mixing the Higgs with a heavier singlet can suppress couplings of 
the h, if it is heavier than the Higgs, the effect is to push the mass of h down, aggravating 
naturalness issues. Thus, we should view a light s as a natural consequence of this model with 
m~ h ~ 110 GeV. Such a scenario is somewhat distinct from previously discussed scenarios 
[41,43] because of the presence of the light s which can be easily produced at rates comparable 
to that of the SM (in general, roughly a factor of 10 smaller). 

Let us further note that we have been discussing the model independent tuning. Whether 
or not one can achieve these mixings in a given model without, e.g., tuning to prevent a 
tachyonic s is a separate, often more stringent constraint. 

As a final comment, we note that these limits are based on the best available limits of 
h — > 2a —>■ 4b, which are still preliminary. Final limits may further constrain this scenario. 

Model Building: Single Stage Cascades into b Quarks 

It is straightforward to construct models in which the Higgs is mixed strongly with a lighter 
scalar. Because of this mixing, if the Higgs decays h — > 2a, then generally 5 decays to 2a as 
well. Here we will discuss the model building and tunings associated with large mixings, and 
situations where h — ► 2a — ► 46. Situations with h — > 2a — > 4g, 4r, where Higgs mass may be 
below 110 GeV, will be deferred to subsequent discussions. 
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Figure 2: In these plots, the contours are of ranges of s mass (a,c) or h mass (b,d). The x-axis 
specifies the value of m s arising from the mo operator. The y-axis is all other contributions to the 
mass squared term for s. The values of ■Jfn^ are 250 GeV for a,b and 350 GeV for c,d. The dashed 
contours are lines of constant sin 2 9 as labeled. The contours (darkest to lightest) are for masses less 
than a) 20, 35, 50, 65, 80, 90, 95; b) 100, 105, 115, 130, 145, 160, 175; c) 20, 35, 50, 65, 80, 90, 95, 
100; d) 108, 110, 115,130, 145, 160, 175. 

There are three operators which can induce significant mixing with the Higgs: mp, A^, 
and Ah. Ah principally mixes s with H° rather than h because cos(a + f3) is small, so we 
focus on the other two. 

rri£) is unique in that while it induces mixing, it also adds a diagonal mass for s, so that a 
tachyon never appears. As a consequence, with this operator, it is simple to get large mixing 
without having to tune masses to a high degree. In figures ^ we see that one can easily achieve 
large Higgs masses and large mixings over broad ranges of the parameter space. 

Xh is somewhat more challenging, because we require a sufficiently large //-term from 
experiment. From searches for the chargino, we typically require [x > 100 GeV [44]. In figure 
H we consider this scenario for A = 0.25. Here we see that it is very challenging to have 
sufficiently large /i while keeping the light state dominantly s (i.e., sin 2 9 < 0.5). One can 
tune this scenario to achieve this, but generally, it is most natural to have sin 2 9 ~ 1 and the 
s decoupled, and m/j > 110 GeV (as is required with no singlet component). Naturalness 
here is not considerably improved from the MSSM, unfortunately. 
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Figure 3: In these plots, the contours are of ranges of s mass (a,c) or h mass (b,d). Here the \h 
operator is taken to be A = 0.25. The x-axis is the value of the effective /x parameter (once s has 
acquired a vev). The y-axis are other soft contributions to the mass squared of s, in addition to 
the X 2 v 2 term arising from A^. The values of y/m^m^ are 250 GeV for a,b and 350 GeV for c,d. 
The dashed contours are lines of constant sin 2 9 as labeled. The contours (darkest to lightest) are for 
masses less than a) 20, 40, 60, 80, 90, 96; b) 100, 115, 130, 145, 160, 175, 190, 205, 220; c) 20, 40, 60, 
80, 90, 105; d) 105, 115, 130, 145, 160, 175, 190, 205, 220. The white region in the lower right corner 
is where s is tachyonic. 



Since A^ is in general more tuned, we list a benchmark with mo below (masses in GeV 
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and in the decoupling limit which we define as via 3> and large tan (3 ^ 5) 
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(5.2) 



Here "tuning" is defined as the tuning in necessary to satisfy the necessary requirements. 
The maximum value of mo arises from the limit on s — > 2a decays and the minimum value 
arises from the lower limit on mi. The tuning of a parameter x which can lie between two 
values a and 6, with b > a will be defined to be the smaller of 100% or (6 — a)/ (b + a). 
Qualitatively, this is the average fractional change in a parameter which is allowed consistent 
with phenomenology and constraints. In future benchmarks, it will generally be defined as 
the largest tuning in <5 2 or <5 2 to achieve the proper spectrum of s, a masses. This tuning is 
usually related to the lighter of the two states. In these scenarios, constraints on Cib+4b ( see 



eq. 5.1) are very stringent, so this region may not be allowed by the full analysis. 



5.2 Higgs Decays Through a Single Stage Cascade (without b's) 

Single stage cascade decays of the Higgs are familiar within supersymmetric theories. Within 
the CP violating MSSM such decays can easily occur [45-47]. Within the NMSSM, the decays 
h — > 2a — > 47 are well studied [36,48], while it has been argued, prior to the combined LEP 
analyses, that h — > 2a — > 46 allowed considerable reduction in the fine tuning of mz within 
the NMSSM [17]. More recently, it has been argued that the combined LEP limits favor a 
h — > 2a — ► 46 decaying Higgs at approximately ~ 105 GeV [41], which we argued in 
section 5A, had applied constraints that may have been too generous, h — > 2a — > 4r has 
also been considered [43,48]. Here we study these, and an additional possibility, namely, 
h — » 2d —> Ag. 

h^2s,2a^4r: Phenomenologically, it makes little difference whether the intermediate 

state here is 2s or 2d, except if it is 2s, the d can be completely decoupled, while if it is 2d, 

the s could mix considerably with the Higgs, giving large corrections to the Higgs couplings. In 

general, models of this sort are highly tuned, because the mass of the intermediate scalar must 

be 1z 10 GeV in order for decays to r's to dominate over decays to b quarks. As a consequence, 

they generally do not lead to a significantly less tuned theory than the MSSM, although they 

offer interesting and distinct phenomenology. 7 Generically, one expects a scalar field with 

mass ^$ 10 GeV, as well as accompanying decays h — * 2s, 2d — > 2t2c and h — > 2s, 2d — > 4c. 

There may also be some interesting B physics in this mass range [50]. 

The only existing analyses on this scenario [18, 19, 22] cut off at > 86 GeV, for 

theoretical reasons related to the MSSM which do not exist in theories with additional singlets. 

The extension of this analysis is important to understanding the true limits on such a scenario. 

h — > 2a — > 4g: One of the most interesting possible signatures results in the 4 glue final 

state. Such a final state is difficult to see at LEP and has not been constrained by any 

7 In our scenarios h decays into 2a are correlated with the size of the a mass which leads to tuning for a to 
be lighter than two b's. If the decays are mediated by derivative couplings, e.g. a Peccei-Quinn axion, then 
this light a can be reasonably natural (see for e.g. [49]). 
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existing analysis, save the model independent analysis [22], allowing such a Higgs, in principle, 
as light as 82 GeV. The extent to which we believe other analyses (e.g., the two-jet flavor 
independent decays) can be interpreted to constrain this scenario will be discussed in the 
benchmark summary. 

Models with this final state can be easily constructed with essentially no fine-tuning, and 
squarks as light as the direct experimental limit. Although a search for four jets at a hadron 
collider seems impossible, the same processes which generate the four jet final state can also 
generate h — > 2d — > 2g2j, which should occur with a rate down by roughly a 2 EM /a 2 ~ 10~ 2 , 
see Mq 1 subsection of section [|. This rate can be increased if couplings of the a to Higgsinos 
exist as well. A search for the two photons at the LHC should be feasible [36], and there may 
be background reduction by requiring two jets which reconstruct to the invariant mass of the 
photon pair (however, this is unlikely at a hadron collider like the LHC). 



Model Building: h -> 2a -> 4r, 4g via X h , Mq 1 



A/j allows decays to both s and a. However, A^ generally mixes the CP even s with h, giving 
strong constraints from s-strahlung. Consequently, with A^ , the most natural models are those 
with h — > 2d — » X. In order to evade constraints on s-strahlung, it is most straightforward 
to simply raise the mass of s through the addition of 5 2 . In doing so, however, one generally 



realizes the large mixing scenario described in section 5.1, requiring that raj > 2ms in order 
to evade constraints. One can raise the mass of s much larger than the other relevant scales, 
as well. In fact, this is quite reasonable, for instance if S appears as the N=2 completion of a 
5-D U(l) gauge multiplet. 5-D gauge invariance protects the mass of a, but not s against tree 
level SUSY breaking terms [51,52]. Thus, for our purposes in this section, we will decouple 
s from the theory. 

For h — > 2d to sufficiently dominate over h — > 26, a suitably large A^ is necessary, and by 
the discussion in section ||, we see that the requirement is A > 0.2. At tree level, this results 
in mass squared of (35 GeV) 2 for 5. 

Let us now consider the subsequent decays of a. If we include the a nonzero Mq , the 
possibility of 5 — ► 2g opens up. However, the presence of Xh requires the presence of Ah, 
which allows 5 — > 2b as well. This operator need not be large, being generated at the loop 
level. At this size 

A h ~ Y^2rn susy ~ .2 GeV (5.3) 

it would be expected to yield a mixing of a and ^4° at the 10 -3 level or smaller, which is easily 
small enough for the decays to glue to dominate. In comparison, in the NMSSM it is difficult 
to have Mq 1 induced decays dominate because the cross term in \kS 2 + XH u Hd\ 2 acts like 
Ah with an S vev and thus fermion decays tend to dominate in this scenario. 



- 22 - 



The benchmark point then has (taking the decoupling limit) 
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Here we have taken fi = 100 GeV, there is no reduction in the effective /i-term since 5 S is 
large. The inclusion of a nonzero Mq 1 as described in subsection Mq 1 of section |] generates 
the decays to two gluons. Due to the kinematics (heavy h, a), there are no constraints on 
these 4 jet decays. 

If Mq 1 = 0, one expects decays to standard model fermions to dominate. For decays of 
taus to dominate requires m„ to be less than the bb threshold. This requires us to tune 5 2 
to cancel the tree- level mass of m 2 ~ (40GeV) 2 . For instance, the above point with <5 2 of 
— (51GeV) 2 achieves the proper spectrum with a 3% tuning. The necessary changes to the 
parameters and the corresponding changes in the spectrum are shown in parentheses in the 
table in the second line, whereas things that are common to the first are not repeated. 

Model Building: h, s — » 2a — * 4g, 4r via mp, A s , Mq 1 

Cascade decays of the CP-even states into 2a can easily dominate with the inclusion of mo and 
A s . The operator A s induces a large trilinear sa? while mn induces 0(1) mixing between s 
and h. Thus if the 2a channel is kinematically accessible, h, s will dominantly decay into it. In 
addition, there is no mixing in the CP-odd sector, thus Mq 1 provides the only decay channel 
for a; hence the final state of the cascade decay is dominantly 4g and, at a subdominant but 
potentially interesting rate, 2g 2j and 47. Notice that in this scenario the vev of s has no 
relation to the fi parameter and the minima are only local minima due to the A s operator. 
As a benchmark point, we take (in the decoupling limit) 
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Some comments on this benchmark are in order. Looking at the first line, it is a large mixing 
point, where s, h are similar in mass and the main constraints are on the h decays. As 
mentioned above, with no Ah operator, there is no mixing in the CP-odd sector so d can only 
decay into gauge bosons through Mq 1 . The constraint from the dominant 4o final state only 
requires m„ ^ 5 GeV [35], so these scenarios are viable and quite natural. 

As a variation on the phenomenology, r final states can also be considered. For the 
4r state to compete, one could include a small Ah operator and tune the 5 mass below bb 
threshold, at the expense of increased tuning. As before the necessary changes in parameters 
are given in parentheses in the second line. The constraints on cascades into r now apply, 
which generally require that mj^ 5 GeV if B(d — > 2r) ~ 1. However, these points are on 
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the edge of a constraint contour, for instance when rag is between 45-76 GeV, one instead 
requires 7 GeV, severely limiting the allowed parameter space [18,19,35]. In summary, 
these benchmark points suggest that cascades into r's seem to be disfavored by both tuning 
considerations and available parameter space, but still remain an interesting phenomenological 
possibility. 

5.3 Higgs Decays Through a Two Stage Cascade 

A very interesting phenomenological possibility is that the Higgs cascade decay proceeds 
through two intermediate scalars. This is a very reasonable possibility in supersymmetric 
theories, because the new singlet actually comes with both the scalar and pseudoscalar, so 
h —* 2s —* 4a — » 4X is possible or with CP mixing, h — > sa — > 3a — » 3X is possible. However, 
when considering theories with A^, it is essentially impossible to achieve two stage cascades as 
the dominant decay mode, see Fig. ||. Furthermore, with A^, s is lighter than 5 due to mixing, 
making it difficult to arrange a scenario where something of the form h — > 2s — > 4a — > 4X 
can occur. 

However, in the broader context we are considering here, such decays can easily arise, 
with only moderate spectral tuning. In general, the models rely upon the presence of D-term 
supersymmetry breaking through the supersoft operator mo- 

There are two forms which we can study, one with CP mixing and one without. With 
CP mixing, we can consider h — > as — > 3a — > 6r, 66. Such scenarios can easily occur with 
essentially no tuning, and occur when the decay h — > 2s is kinematically forbidden. In cases 
where the six b-quark final state dominates, we also expect a non-trivial fraction of 46 2r 
events as well. 

Without CP mixing, the cascade generally follows h — > 2s — > 4a — > 8g, 8r, 86. Scenarios 
with so many soft final states have not been studied at LEP, and thus have no constraints 
beyond the OPAL model independent constraints (nih > 82 GeV). This allows very light 
stops, and the only tuning required is in achieving the proper spectrum of scalars. 

In cases with eight gluon jets, there are typically levels of 6527 which are non- negligible. 
However, the photons are sufficiently soft that the background is very large, and the jets are 
very soft, well below the cuts to be applied at the LHC. Likewise, the decays to eight b quarks 
or eight r's seem difficult to study at hadron colliders. These scenarios probably will instead 
be studied by a future linear collider like the ILC. 

Model Building: h -> 2s -> 4a 8g, 8r, 8b via m D , A s , Mq 1 

With the operators rri£i,A s , and Mq , it is possible to find points that have significant two 
stage cascades. With the following benchmark point, the Higgs-like mass eigenstate is pushed 
heavy by mixing with the s, and via dominance of mo over A s , decays h — > 2s are favored 
over decays into 2a. The small value of A s is still sufficient to have s decay into 2a's which 
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completes the two stage cascade. 
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At this point, the constraints that are independent of the a decay products are: 1) decay 
independent analysis on s which requires m% > 5GeV for the value of sin 2 9, 2) SM searches on 
the S and h which are both satisfied with stronger constraints on h since £| M = cos 2 9B{h — > 
SM) = .14. 

The final state products of a determine what additional constraints apply to this scenario. 
In the first line of the benchmark, we will consider gluon and b quark final states. With 
turned off, a nonzero Mq 1 causes it to decay primarily to gluon jets. For this final state, 
there are no further analyses that constrain this scenario. However, if the a does mix with 
A by both turning on a small Ah operator and reducing Mq 1 , decays 6's dominate. In this 
case, the limits apply to the effective £f 6 + 46 defined earlier. We find that these are consistent 
with the new LEP limits, assuming that the analyses are not sensitive to the 86 final state 
(this sensitivity has not yet been determined by any LEP analysis). If a is lighter than 10 
GeV, again listed in the second line with new parameters/changes in parentheses, r's can 
dominate, and is consistent with the new LEP analysis of 4r events from s-strahlung and 
Higgs-strahlung. Finally, in the case that Mq 1 and mixing decays are comparable, more 
complicated scenarios are allowed, with many possible final state topologies. 



Model Building: h — > as — ► 3a — ► 6b, 6t via rnfyp, tud, A s 

In the presence of m/Qp, and A s , the coefficient of the hsa coupling is 
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Where we have taken the decoupling and large tan (3 limits, as well as the mj = limit and 
at tree- level A = gy + g 2 ■ 

We are most interested in scenarios where 2ma < m§ ~ 60 GeV. Using the perturbative 
expression for sin 9^ (eq. |4.27j ), we require sin 2 9 a h < 0.015, which is a typical upper limit for 
a — * 2b. The constraints on a — > 2r are roughly sin 2 9 a h < -07 in the kinematically allowed 



region. Comparing this limit with eq. 4.27, we have (neglecting terms down by m~/m~) 

m 2 CP < .12 mj/ (sin 9 cos 9) (5.8) 
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for b quarks, and 

m 2 CP < .26 mf /(sin 9 cos 9) (5.9) 

for tau decays. Hence, m^p can be 0(m~), and thus not a small parameter. One can still use 
the approximate expressions for estimates and intuitive understanding, but for better than 
0(1) precision, we must calculate exactly. 

The presence of A s is important, both for h — > as and s — > 2d decays. In general, one 
does not need a large value for A s , with 4 GeV sufficient in order to secure sufficiently large 
h — > as and s — ► 2d decays. 

Let us consider the benchmark point (mass units are in GeV): 
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In line one, we consider b quark decays (gluon decays cannot dominate since there is mixing 
between a and h). This point is consistent with the limits on the effective £|&+4b f rom 46 and 
2b decays. To get to r decays, one changes the parameters as listed in line two, which satisfies 
the constraints on 4r rates. Again, these results assume limited sensitivity at LEP to the 6 
final state decays (6b and 6r). 

6. Phenomenology and Benchmark Point Summary 

It is remarkable that by simply extending the MSSM to include a singlet superfield, there 
can be so many modifications to the Higgs phenomenology. Furthermore, many of these 
modifications do not arise in the NMSSM framework. 

To aid in understanding, we have attempted to isolate a few points in parameter space 
which demonstrate the relevant phenomenology. There are several qualitative issues which 
are relevant. Let us begin by reviewing our different decay scenarios. 

i) h — > 2a — » 4b: The new LEP limits are noticeably more constraining, and require a certain 
degree of suppression in the Higgs production and branching ratio of a to b quarks. This can 
be achieved with ^/rn^~JTT]~ = 300GeV (O(10 — 15%) tuning). In a crude attempt to combine 
limits on 2b and 46 decays, we defined a cumulative t^b+Ab ( m ec l- H3) an ^ required it to be 
less than y/2. This is allowed principally in the "just so" region, where we have essentially 
saturated the LEP limits, i.e., pushed the parameters so that LEP was nearly sensitive to this 
scenario. Hence, these regions may be impacted by the finalized LEP limits. Hidden tunings 
seem to appear when we attempt to use Aft,, the operator used in the NMSSM, to generate 
the decays because of the large [i term and the subsequent large mixing term with s. Such 
hidden tunings are smaller in models with the supersoft operator rri£, together with A s . 

ii) h — > 2s/2a — > 4r: We have not demonstrated a model which can achieve these decays 
above the kinematical threshold (i.e., where a — > bb is not kinematically forbidden). However, 
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we can achieve these decays by tuning the models to the 10 percent level (necessary to get 
the a mass into this region). Furthermore, it is troubling that the best available limits in 
this region (from OPAL [35]) seem to stop at 86 GeV, a theoretical prior due to constraints 
in two Higgs doublet models which does not apply in cases with singlets. This motivates a 
reanalysis of the LEP data without the theoretical bias that m h < 86 GeV if ma ^ 10 GeV. 

iii) h —* 2a — > 4g: It is remarkably simple for this decay to dominate under the assumption 
that Ah is loop suppressed or absent. Such models could be very natural with arbitrarily 
light stops (subject to direct search limits), and a Higgs with mass m h as light as 82 GeV, the 
limit from the OPAL model independent analysis [22]. These scenarios can arise easily with 
either or A/,, and generally come with associated h — ► 2a — > 2g 27 decays at the 7 x 10 -3 
level. 

iv) h — > as —* 3a —* 6b, 6t: Such a decay can dominate Higgs decays via the inclusion of the 
rriD operator, incorporating the sensitivity of the Higgs to D-term supersymmetry breaking, 
but is difficult to engineer with only A^ (i.e., without D-term breaking). This is an important 
example of phenomenology which would not occur in the NMSSM, but easily arises within a 
general operator analysis. Here there is no analysis to exclude this scenario, meaning a Higgs 
as light as 82 GeV is allowed. In the case that the final state is composed of b-jets, such a 
scenario can be quite natural. In the case where the final state is composed of r's, the tuning 
is significant, again at the few percent level. Such a scenario would be difficult to detect at 
the LHC. Finally, although this scenario seems CP-violating, it can be consistent for a, § to 
be both CP-even and thus there are no additional contributions to CP-violating observables 
such as edms. 

v) h —* 2s — ► 4a —* 8b, 8t, 8g: This scenario, too, can only dominate with D-term breaking, 
and not within the narrow NMSSM framework. With D-term breaking, it appears to be 
necessary to have roughly 10% tuning in order to achieve the proper spectrum (with b or 
g final states - a few percent if the final state is r), but again a light Higgs (as light as 82 
GeV) is allowed (however, in our benchmark, an additional 46 rate pushes up the required 
Higgs mass). Because the final state particles are so soft, it is difficult to envision a scenario 
in which the LHC could detect this Higgs. 

The benchmarks realizing these phenomenologies are summarized in Table ^. The bench- 
mark points illustrate the importance of the general operator analysis. Some scenarios are 
only natural with the presence of D-term SUSY breaking. The gluon final states only occur 
once we additionally consider the effects of new, heavy fermions. Such effects are typically 
excluded from NMSSM analyses and yet we see they can generate some of the most interesting 
phenomenology. Despite this, one should keep in mind that unexplored operators may also 
generate the same phenomenology, and thus the phenomenology presented should not only 
be considered in the context of the particular model realizations. 
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Table 2: Benchmark Point Summary. Terms in parentheses are variations on the spectrum in order 
to achieve a — ► tt decays. When A/j is labeled 0/small, the indicates the presence of nonzero Mq 1 
allowing a — > gg decays, while the "small" is to indicate a small level allowing a — > 6fe or a — > tt 
decays, without modifying the spectrum, [as] indicates that the decay is h as rather than h — > 2s. 

6.1 Possible exclusions in existing data 

To the extent that many of these signals have not been explicitly analyzed, one can argue 
that only the model independent bound from OPAL truly limits them. However, once one 
specifies a given decay mode, it is almost certain that the bounds will improve from the model 
independent limit. 

Going further, extrapolating from existing limits might give us an estimate on the poten- 
tial exclusions of an actual analysis. For instance, in searches for h — > 46 in events where the Z 
decays hadronically, DELPHI and OPAL force the whole event into 4 jets, which is reasonably 
efficient even though the total process can have more jets in it (up to 6). So since this anal- 
ysis is akin to the SM 4 jet analysis done for h — ► 26, Z — > hadronic, the flavor-independent 
analyses may be used to estimate the potential limits on h — > Ag. If one assumes that the 
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efficiencies to reconstruct the 4g and 46 state are the same, and that the background for 4g 
is the same as for 2g, one can estimate the scale of exclusion that may be possible in existing 
data. Doing this, it appears that h — ► 4j below 86 GeV has strong exclusion, that interesting 
constraints could possibly be set for 86 GeV < < 100 GeV, and for m~ h > lOOGeV it 
appears unlikely that interesting constraints could be set. All of this is an extremely rough 
estimate, but suggests additional analyses in the hadronic channels should be done. 

Decays of the Higgs to six and eight parton final states are more difficult to estimate, for 
instance it is not known how sensitive the current 26 and 46 searches are to actual 66 or 86 
decays. Analyses that attempt to utilize the multi-jet nature of the decays might be useful, 
but might encounter issues like jet-finding algorithms constructing fake jets in the signal or 
background (OPAL and DELPHI use the DURHAM algorithm which is known to have such 
issues). It is also not clear how much b-tagging can help in the cases of multi-b decays. 

6.2 Future experiments 

It is not clear to what extent these Higgses may evade detection at the LHC. The six and eight 
parton decays almost certainly will be challenging. Similarly, the decay of the Higgs into four 
gluon jets is a tremendous challenge at any hadron collider. The decay h — > 2a — ► 2g2j should 
be more promising. Especially as loops involving the Higgsinos can amplify this beyond the 
expected 7 x 10~ 3 branching ratio, analyses of such scenarios, particularly at the Tevatron, 
would be highly motivated. The associated h — > 4j decay has the added signal of seeing 
both the two a's and the Higgs that decayed into them. Unfortunately, the expected 10~ 5 
branching ratio is probably too small to be seen, so this channel requires either an enhanced 
rate either through increased branching ratio or Higgs production via colored sparticle decay 
chains. 

7. Summary and Conclusions 

One of the strongest predictions of the MSSM is the character of the Higgs boson, both in its 
mass and its couplings. However, the simplest extension of the MSSM, adding a singlet, can 
dramatically alter the Higgs boson phenomenology, in particular by introducing a remarkable 
set of decay scenarios. We have studied, in a relatively general sense, the effects of singlets 
on the decays of the Higgs boson in supersymmetric theories. There are a number of general 
points to be gleaned from this analysis. 

First of all, the MSSM expectations of the Higgs boson decays can easily be modified. 
The presence of the singlet opens up a number of cascade decay possibilities, which are much 
harder to constrain. Secondly, the NMSSM (where the singlet acquires an expectation value 
to generate the //-term) is too narrow a framework to realize many of the interesting decay 
scenarios. A number of them can only dominate the Higgs decay once one includes the 
possible effects of D-term vevs on the physics of singlets. Some cascades can occur but are 
remarkably tuned within the NMSSM, but are not tuned in a broader framework. This makes 
it essential that we continue to include these operators in future analyses. 
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Finally, while these models admit much lighter Higgses, and hence much lighter stops, 
than the MSSM, the sensitivity of mz is typically inadequate as a measure of tuning. Often 
the most severe tuning comes from finding an appropriate spectrum such that cascades can 
occur, while still maintaining a large coupling to the Higgs boson. All analyses of cascade 
decays must be careful to consider these additional tunings rather than simply focus on the 
tuning of mz- Also, through mixing, heavier Higgses can be allowed for lighter top squarks 
and hence with less tuning of the mz- 

Such issues notwithstanding, it is quite straightforward to achieve models which are not 
tuned, and have a Higgs below the MSSM LEP limit. Cascade decays h — * 2d — * 46 allow 
lighter, more natural Higgses, but at the expense of typically saturating the LEP bounds, a 
tuning of sorts on its own. Plus, these allowed regions will likely be impacted by the finalized 
LEP limits. Other cascade decays, in particular h — ► 2d — > 4g and h — > as — ► 3d — ► 66 can 
occur with very little tuning of the d mass, and light Higgses (mr > 82 GeV). Decays with r 
final states {h — > 2d — > 4r, h — > as — ► 3a — > 6r, h — > 2s — > 4a — > 8r) seem to be quite tuned, 
since there is no apparent symmetry to explain the very light a (except in case of PQ axion, 
see footnote ^). On the other hand, two stage cascades h — > 2s — ► 4a — > 86, 8g can occur with 
moderate tuning. In our realizations, two stage cascades (those with six- and eight parton 
final states) require D-term breaking while gluon final states arise with new, heavy fields, 
although there are possibly additional means to produce this phenomenology. 

The collider tests of these scenarios are largely unexplored. Firstly, since the model inde- 
pendent decay analysis of OPAL [22] is the dominant constraint for most of our scenarios, it 
would be extremely useful for a LEP-wide analysis to be carried out. LEP may also have been 
sensitive to some of the new non-standard decays we have introduced, which motivates further 
investigation by the LEP Higgs search collaborations on these specific decays. Therefore, we 
strongly suggest that the following analyses be performed: 

• h — > anything (model independent). Such an exclusion performed by the OPAL analysis 
is an important consideration for all models of non-standard Higgs decays, and thus 
should be performed by as many LEP collaborations as possible. 

• h — * hadronic (i.e., 4 jet, or multi-parton final states). Given that photonic, leptonic, 
and two jet decays of the Higgs are already highly constrained, it is worthwhile to close 
a significant means by which a Higgs might hide, namely in a generic hadronic decay. 
However, such an analysis may be quite difficult to separate from background, even in 
events where large amounts of 6-flavor is required. The simplest analysis should be an 
extension of flavor-independent Higgs search to flavor-independent cascades. 

• h — > 2a — > 4r above = 85 GeV. The theoretical considerations which appear to 
have truncated the analysis at this mass range do not apply in general models, and 
likely significant constraints can still be placed above this value. 

At the LHC and the Tevatron, the most promising decay channel is the h — > 2a — > 2g 2j 
that is associated with 4 gluon decay. At the Tevatron, such a search may be possible, because 
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the jet threshold is sufficiently low. At the LHC, in this topology, the a photon decay could 
be seen, although a more careful analysis is needed [36]. Since the two gluon jets would be 
difficult to measure/choose, this topology does not allow detection of the Higgs. The 4 photon 
decay is probably at too small of a rate to detect, although this may be increased either by 
increasing the branching ratio or if Higgses appear in cascade decays of colored superpartners 
(and thus boosting the overall production rate). Many of the other channels seem even more 
difficult to search at the LHC, but on a positive note, at the ILC, one should be able to detect 
the Higgs through the Z recoil method (see for example [53] and for a linear collider analysis 
of a SUSY model with nonstandard Higgs decays see [54]). 

It should be noted that while many of these scenarios are very difficult in terms of Higgs 
detection, they are not nightmare scenarios. Quite the contrary, the existence of these decays 
allows a light Higgs and would naturally be associated with a wealth of light superpartners. 
This suggests that the proper lesson of LEP is not a need for models with a heavy Higgs, but 
rather that further thinking about the possibilities of a stealthy Higgs is in order. 
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A. Trilinear couplings 

rriD — Supersoft Operator 

The terms arising from the supersoft operator are 
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while those arising from the D-terms are 
D-terms 
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where at tree- level A = §y + g 2 ■ Where the notation is as in section ||; is the mixing angle 
between s and h and (p is the mixing angle between a and ^4°, from the operator. 

Top corrections 

The loops from the top induce trilinear terms with an overall coefficient C = log (m^m^/ra 
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The Xh and fi term give trilinears 
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In this case the induced trilinears are 
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This soft term gives the following trilinears 
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Since we are ultimately interested in h — > sa and s — ► aa decays, we find those terms in the 
Lagrangian. Arising from the usual MSSM D-terms, we have (at leading order in m 2 CP ), 
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While from mo we have 
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This term is often important in conjunction with A s . The equivalent trilinears with the 



combination of A s and rr? cp are 
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